Optical injection and terahertz detection of the macroscopic Berry curvature 
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We propose an experimental scheme to probe the Berry curvature of solids. Our method is 
sensitive to arbitrary regions of the Brillouin zone, and employs only basic optical and terahertz 
techniques to yield a background free signal. Using semiconductor quantum wells as a prototypical 
system, we discuss how to inject Berry curvature macroscopically, and probe it in a way that provides 
information about the underlying microscopic Berry curvature. 



Berry's phase permeates many fields of physics. In 
quantum mechanics, the net phase acquired by a wave- 
function during a cyclic change in the parameters of a 
Hamiltonian is a gauge invariant quantity that has mea- 
surable effects. This quantity can be directly expressed 
in terms of a physical property, called the Berry curva- 
ture (BC), which relates to the Berry phase in much the 
same way as the gauge- invariant magnetic field relates 
to a gauge-dependent vector potential [lj. In materials 
science, it is rapidly becoming clear that the role of the 
BC is as fundamental as that of an energy band [2 J . The 
representation of position operator in terms of the Bloch 
states of a crystal is intimately related to the BC, and 
it thus appeared in the works of Adams and Blount |3J 
on the topic much earlier than the rigorous formulation 
by Berry |4J. In 1990s, King-Smith and Vanderbilt intro- 
duced a theory of electric polarization of solids as a bulk 
quantity [5j, the present version of which is grounded 
in Berry's phase [2, 6J. In parallel to this, the BC has 
appeared as a central physical quantity in the work on 
anomalous Hall effect (AHE) in ferromagnetic materials 
[7J [8], and the intrinsic mechanism of spin Hall effect 
(SHE) [9} [TO]. Haldane expressed the non-quantized part 
of the intrinsic Hall conductivity in terms of integral of 
the Berry connection (vector potential of BC) on the 
Fermi surface [llj, thus recapturing the essence of Lan- 
dau's Fermi liquid theory. Recent works also show the 
role of BC in the photogalvanic effect p~2j [13]. 

Central to these effects is anomalous velocity, which 
refers to motion of charges perpendicular to their usual 
group velocity |3J. As an average of the BC it appears as 
its main experimental manifestation. Though the above 
works elucidate the role of the BC in explaining various 
phenomena, they are severely restricted as methods to 
probe it as a basic property of solids. This is because 
the d.c. response reduces to an integral of Berry cur- 
vature over the full BZ, and for partially filled band, to 
the volume occupied by the carrier distribution function 
[TT] , The measured response is thus restricted either to 
the full BZ, or slightly displaced equilibrium occupation 
function of a partially filled band. 



This letter presents a method to study this important 
and complementary quantity of solids, which is much 
less restricted in its sampling of the Brillouin zone. It 
was anticipated by the work of Moore and Orenstein[T4], 
who showed how confinement in quantum wells induces 
a Berry's phase that leads to a helicity dependent pho- 
tocurrent. While we focus on the hole bands of a GaAs 
quantum well (see Fig. [T]), chosen because of their large 
BC |9J with a rich structure, our method is very general 
and in principle can be applied to a wide range of mate- 
rials. An optical excitation by circularly polarized light 
induces a non-zero transient macroscopic Berry curva- 
ture (TMBC) in such materials, due to the creation of a 
state that breaks time-reversal and space inversion sym- 
metry. The injection is followed by a linearly polarized 
terahertz wave, which drives the optically injected carri- 
ers. In the presence of TMBC, the motion of charges has 
an anomalous component, which is perpendicular to the 
polarization of the terahertz wave. The detection of the 
radiated THz field in the cross-polarized direction is thus 
entirely a signal due to the TMBC. Such a TMBC can be 
injected even in materials of Chern class zero, where an 
equilibrium MBC is not allowed [15J. This is in contrast 
to spin hall insulators in which a net MBC exists in the 
ground state. 

Our letter is organized as follows. We begin with a dis- 
cussion of the microscopic Berry curvature. We then dis- 
cuss the optical injection of TMBC, its detection by lin- 
early polarized THz excitation, and its intrinsic lifetime 
due to electron-hole scattering. Finally, we present the 
results of our numerical calculations of the microscopic 
and optically injected TMBC, followed by the results of 
the anomalous velocity induced by the THz excitation. 



MICROSCOPIC BERRY CURVATURE 

Bulk crystals: The band theory of solids de- 
scribes single electron wavefunctions by Bloch functions, 
i(j n (k;r), for energy bands, e n (k). Here n labels a single 
energy band with a continuous slope, starting from the 
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and cr 3 is the third Pauli matrix. Figure [2] shows the top 
two valence energy bands and the corresponding Berry 
curvature for a 15 nm thick, quantum well. The middle 
panel shows Q^ n (k) defined in ([3|; large Q^ n (k) results 
from a large mixing of two or more states [2 lj, here arising 
because of contributions from the light hole to the two 
bands. 



MACROSCOPIC BERRY CURVATURE 



Figure 1: Illustration of the proposed scheme for injection and 
detection of macroscopic Berry curvature 

energy levels at the center of the Brillouin zone [16J. Be- 
tween two degenerate bands, or within a single band, the 
vector 

Cnm( k ) =i J dru* n (k;r)-^u m (k;r), (1) 

acts like a vector potential in the momentum space dy- 
namics of electrons, its dependence on the choice of the 
phases of the {u m (k;r)} leading to the analogue of a 
gauge dependence. Taking the curl of Cnm(k) l ea ds to 
the (microscopic) BC, a "gauge invariant" quantity at 
each k, 

Km(k) = (Vx(nm(k)r-ie a01 £ [<£,(*), Q m (k)] , 

(2) 

In the presence of time-reversal and spatial inversion 
symmetries, Cl(k) takes the form of a traceless matrix 
within each degenerate subspace |3J; the macroscopic 
Berry curvature (MBC), 

(ft) ^ ttnm(k)pmn(k), 
nmk 

where p nm (k) is the single particle density matrix, then 
vanishes in equilibrium; even in a material such as GaAs, 
where there is no inversion symmetry, the equilibrium 
MBC vanishes because the band structure is of Chern 
class zero. Away from equilibrium, such as in the d.c. re- 
sponse of a p-doped semiconductor, the MBC can play a 
role [8-10J, but optically excited distributions allow much 
greater access to its local probing in k-space. 

Quantum wells: We focus on the valence states of 
a quantum well, described by the Luttinger model of a 
square quantum well grown along the [001] = z direction. 
In this model [17] there are degenerate (±) wavefunctions 
described by 4-component spinors, labeled f^(k;z) for 
each 2-dimensional subspace n, and the microscopic BC 
is equal to z^ n (/c)<7 3 , where 

Kn(k) = e zia , f dzd^f-\k;z)-d u f-(k;z), (3) 



No macroscopic effect results from Q^ n (k) unless an 
imbalance is created between the hole populations of 
the / ± states of these subspaces, leading to a nonzero 
(Q z ). In optical excitation across the bandgap with left- 
circular ly polarized light, the matrix elements connect 
/~ to the spin up conduction subbands with a prob- 
ability 1/3 of the corresponding transition from to 
the same conduction subband. This yields a 3:1 popu- 
lation imbalance between the two otherwise degenerate 
states, leading to a non- vanishing (Q z ). The possibility 
of creating this imbalance in the non-equilibrium state 
is a well-established experimental fact [18J. A correction 
arises from the electron-hole interaction such that the 
final state is an exciton. This correction is small for ion- 
ization states (the only states relevant above the optical 
gap), and decreases the higher up we move in the band. 

To access (fl z ) we imagine a linearly polarized THz 
electric field, E(i)i, driving the optically injected hole 
population. Terahertz excitation does not couple and 
/~ states since there is no momentum matrix element be- 
tween them, and therefore it does not disturb the injected 
(Q z ). The Schrodinger equation for the hole density ma- 
trix p(k, t) coupled to the external THz field polarized 
along x (or [100]) direction reads, 

= [H (k) + eE*(t)C(k), p(k, t)] + p\ scatt . 

Here Ho(k) is the Luttinger Hamiltonian, and the last 
term is the contribution of Coulomb scattering. A gauge 
transformation exists such that ( x (k) = locally in the 
region of interest, and within this gauge we can easily 
compute the solution of Q with p(fc, 0) set equal to the 
hole populations created by optical excitation. The ex- 
pectation value of the velocity operator, v, is then calcu- 
lated at each time point to study the macroscopic motion 
of charges, 

(v(t)) = ^2 Pnm(k, t)v mn (k). (5) 

k 

Sum rules that relate ft to £ [16J can be used to show 
that in the linear regime an anomalous contribution to 
the average velocity may exist, and is proportional to the 
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Figure 2: (a) Highest valence (heavy) hole bands for a 15 nm thick [001] GaAs quantum well; (b) Berry curvature Q z (k) (lines 
as in (a)); (c) (Q z ) (thick solid line) and (O^) for odd parity states (lines as in (a)); (d) Coulomb scattering rate between the 
+ and — states for optical excitation at 90 meV above gap energy Lines as in (a) 



macroscopic Berry curvature, 

(v(t)) anomalous = -yl(W)E*(t), (6) 

until (Q z ) decays due to scattering between the ± state 
populations, primarily through the dynamic polarizabil- 
ity of the electron gas, described by the term p\ scatr The 
electrons generated by the optical excitation would al- 
ways diminish the Berry phase effects, and therefore this 
is an intrinsic lifetime because it survives in the limit of 
no impurities. The self energy effects may be decomposed 
into two contributions: (A) decoherence within each set 
(either + or — ) of states, and (B) population transfer 
between the sets of states. Type A limit the sensitiv- 
ity of the dynamics within each set to the microscopic 
BC, while type B reduce the net effect of MBC; B ef- 
fects degrade the signal more effectively than A effects 
by restoring the time-reversal symmetry broken by the 
initial optical excitation. 

We model these effects by setting pnm(k,t)\ scatt = 
—T] nm p nrn (t) for off-diagonal terms (type A; n 
and m in the same set), and Pnn(k,t)\ scatt = 
-7nn(fe) [Pnn(k, t) - p n n(k, i)] for population relax- 
ation terms (type B; here n(n) refers to the f^(fn) 
state). To estimate the rates j nm (k) and 7 n n(fe), we 
start with the RPA self-energy, employ the generalized 
Kadanoff Baym ansatz |19| , and follow by the Markov 
approximation. The Auger scattering of holes by 
electrons is negligible, due to the energy conservation 
condition at each vertex in the self energy diagram. 

RESULTS AND DISCUSSION 

For excitation with a 100 fs Gaussian optical pulse, the 
resulting (Q z ) is shown as a function of energy above the 
band gap in the third panel of Fig. [2j The panel also 
shows the breakdown of (Q z ) into contributions from the 
individual bands. The initial rise in (Q z ) is entirely due 
the first band, in which the dominance of / + component 
leads to the sign opposite to its microscopic counterpart 



in the middle panel. The large contribution from the 
second band starts to dominate above 85 meV. 

The rightmost panel in figure [2] shows the scattering 
rate 7 n ,n(fc) f° r eacn of the two bands at a carrier den- 
sity of 10 11 cm -2 . The rate remains below 1 THz, and is 
suppressed at small where the spinors have orthogonal 
dominant contributions. We find that the rj nrn (k) range 
between 1 to 4 THz, but their effect on the THz dynam- 
ics is only secondary as discussed in relation to type A 
interactions. So we can expect that (Q z ) does not vanish 
during the subsequent THz probe. Thus we see that the 
effects of relevant scattering between electrons and holes 
are significantly suppressed by symmetries of the well. 

We now turn to the results of our numerical solution 
of Q, and present the normal and anomalous hole ve- 
locities calculated via The field is polarized along 
x so that the anomalous velocity is along y; see Fig- 
ure 1. Our calculations are performed with E(t) = 
E e-^- to)2/2r p cos(uj t). We chose E = 0.1 kV/cm, 
uq = 1 THz, r p = 1 ps, and to a conveniently chosen point 
in time. By comparing the velocities with and without 
the p\ scatt term in Q, we found that the the anomalous 
velocity is reduced by c± 36% while the normal velocity 
undergoes little change by scattering as expected. 

In Figure [3] we demonstrate the relationship between 
the magnitude of (v y (t)), and the hole populations 

for different photon energies above the gap; large anoma- 
lous velocities clearly arise when the Berry curvature of 
the populated states is large. These calculations lie in 
the linear regime, where ^ is a good approximation; in 
the non-linear regime, no clear cut relationship exists be- 
tween the expectation value of the velocity operator and 
the (Cl z ) as the concept of intraband motion itself breaks 
down. 

As shown in Fig. [T] the anomalous velocity could be 
detected by measuring the emitted THz radiation polar- 
ized perpendicular to the incident field. While the much 
lighter electrons in the conduction band would also emit 
THz radiation, and at much higher power than the holes, 
the Berry curvature of conduction bands in these quan- 
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this macroscopic effect is at least a few picoseconds, and 
its robustness results from the properties of hole wave- 
functions under the symmetry operations of the quantum 
well. We have presented a scheme to make this macro- 
scopic effect accessible experimentally via the anomalous 
contribution to the THz emission from a quantum well 
driven by a linearly polarized THz field. A successful im- 
plementation of this scheme would open up a new venue 
in exploration of the Berry curvature as a fundamental 
property of solids. 

The authors acknowledge financial support from Natu- 
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Figure 3: THz induced velocities (left) and populations (right) 
at 20 meV, 90 meV, and 350 meV (top-bottom) for 15 nm 
wide QW. Anomalous (normal) velocities are in solid (dashed- 
dotted) lines. 



turn wells vanishes. So only the anomalous velocity of 
holes would contribute to the THz emission perpendic- 
ular to the incident field. In addition, placing the THz 
detector at an oblique angle in the plane of the incident 
field vector, the dipole field of the parallel component 
can be highly suppressed while leaving the perpendicular 
component unaffected. Separation of the desired signal 
from background is thus already built into this method. 

In the above analysis, we only took intrinsic scatter- 
ing due to Coulomb interaction among electrons into ac- 
count. We emphasize that a large hole mobility is cru- 
cial to the experimental success of this scheme. Spin- 
flip scattering of holes is an approximate concept in this 
scenario, but it is a useful characterization of the detri- 
mental effects of impurity scatterers in our scheme. It is 
suppressed exponentially close to the T point, but rises 
sharply to rates faster than 1 THz as a function of k 
for areal concentrations of 10 10 cm -2 [20J. The effects 
of phonons are subtle. On the one hand, they may pro- 
vide a dominant scattering channel. On the other hand, 
as observed previously, only those scattering events that 
link the partner states are actually detrimental. Detailed 
calculations of these effects, including lattice vibrations, 
is the focus of our ongoing work. 

In conclusion, we have shown that significant macro- 
scopic Berry curvature can be injected in GaAs quan- 
tum wells by circularly polarized light. The lifetime of 
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